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We discuss the various physics aspects of hypothetical, coherent excitations of the 
pion field caused by spontaneously generated local distortions of the chiral order 
parameter. Such distortions, which may occur in the expansion of a hot fireball 
created in relativistic nuclear collisions, may be influenced or even triggered by 
the interaction with strong electromagnetic fields through the axial anomaly. We 
suggest that semiperipheral collisions may provide the best environment for the 
formation of these disoriented chiral condensates. 



1 Introduction 

Pion emission with an asymmetric charge ratio, N(ir°) ^ | (N(tt + ) + N(w~ 
may result from the decay of a coherent, isosinglet excitationxifthe pion field 
Such an excitation can be formally represented in the fornmu 
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Here a\ denotes the isovector creation operator for a pion with momentum 
k. Specific charge states are indicated by a superscript (0, ±) and <p(k) is the 
coherent wavefunction in momentum space. A simple geometrical argument in 
isospin space U shows that the ratio of neutral to all pions contained in such a 
state follows the probability distribution 

P{R) = -±= where R = N(ir°)/N(ir). (2) 

while the average charge ratio (i?) = f Q P{R)RdR has the expected value ^, 
the distribution P(R) differs strongly from the Gaussian 

P stat (i?)~exp(-fc-^) (3) 
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that holds for a statistical ensemble of pions with net isospin zero. Essentially 
the same result is obtained for a coherent 1 = state containing 2N pions, for 
which the probability of finding 2n neutral pions is 



n\ ) (2JV + 1)! 

. One 

imitO 



in the limit n, N — > oo,R = n/N fixed. One also finds strong correlations 
among the charge channels; in the same limit I 

4 2 1 

Coo -> r, C + --, C++ -> -, etc., (5) 

where 

(jV(7r Q )jV(7r 6 )) 
~~ {N(w a )}{N(7rb)) — 1 W 
Clearly, neutral and charged pions are anticorrelated, whereas pions of equal 
charge (especially neutral pions) are strongly correlated. 

2 Formation of DCC 

What makes these formal results interesting in the observationi'i that such co- 
herent, isosinglet pion excitations could occur spontaneously in the de-excitation 
of a highly excited region of space, in which the chiral order parameter - the 
quark condensate (qq) - has been restored to its symmetric value (qq) = 0. As 
the order parameter returns to its equilibrium value | (qq) \ — f„ , its orientation 
in isospin space may temporarily differ from the isoscalar, true equilibrium 
configuration 

—=(uu + dd) = /„■, (ud) = (du) = —={uu — dd) = 0. (7) 
V 2 V 2 

Hence the name disoriented chiral condensate (DCC). The formation of a lo- 
cally homogeneous, albeit disorientated, region of chiral condensate is analo- 
gous to the growth of a domain of local magentization when a ferromagnet is 
quickly cooled ("quenched") below the Curie temperature. 

This process can be conveniently analyzed in the framework of the linear 
sigma modelQ Denoting the chiral order parameter field by ((f>o,(f>) 

4>o = (qq), (f> = (0q), (8) 
one finds that a mode of momentum k follows the equation 

^$ k (t) = (\fl k* - A(0(i) 2 )) (9) 
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The order parameter grows exponentially when the factor in parentheses on 
the right-hand side is positive, i.e. k 2 < A(/ 2 — {4> 2 ))- This condition requires 
that 

1. the mode has a long wavelength and 

2. that the expectation value and the fluctuations in the order parameter <j) 
are small compared with /„■. This condition is fulfilled after a quench. 

Although the exponential growth affects modes with small k, it is not easy 
to grow spatially large domains, because the growth shuts itself off when (4> 2 ) 
approaches the equilibrium value f v . One finds that domains are characteris- 
tically of size r = ("v/A/tt) -1 ~ 0.5 fm for A = 20 demanded by sigma model 
phenomenology!!! It turns out that it is more favorable to consider a region 
of space in which some thermal fluctuations of the order parameter remain. 
Denoting the effective potential as 

V eS (ct>) = -\a 2 ^ 2 + \\4>\ (10) 

where /i 2 is a function of T, with /z 2 (0) = A/ 2 and fi 2 {T c ) = 0, one finds that 
the domain growth is governed by the correlation function u 

C(r, t) = (0(0, t)<t>(x, t)) exp (at - (II) 

where r = |x|. As T approaches T c , n becomes small, implying slower expo- 
nential growth, but also longer correlation lengths. As growth stops after a 
time tf ~ A* -1 , the domain size is (r 2 ) 1 / 2 ~ H~ x . 

In the real world, since chiral symmetry is explicitly broken - correspond- 
ing to the presence of an external magnetic field in the magnetic analog - no 
disorientated domains will form when the system is annealed. Very small do- 
mains will form if the system is fully quenched. The ideal scenario must be 
somewhere in between. It turns out that a three-dimensional expansion of a 
hot region of space filled with pion gas is close to the ideal conditions lij As the 
region expands and cools below T c , the long wavelength pionic modes enter 
the unstable region but the system stays close to the critical line n = 0. This 
expectation is confirmed by numerical calculationsEH 

3 Decay of DCC 

The coherent pion excitation eventually decays into individual pions by deco- 
herence of its individual components of different momentum due to residual 
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interactionsE3 Consider a coherent state in the 0(4) linear sigma model with 
explicit symmetry breaking 




(12) 



which has "rolled down" to a minimum of the potential along the chiral circle 
\4>\ w / w . Neglecting explicit symmetry breaking for a moment, the radial and 
angular directions at this held configuration define isoscalar (a) and isovector 
(71") modes in the local environment of the DCC. The decoherence of the co- 
herent state is facilitated by the special form of the interaction between the 
radial and angular modes of oscillation. Consider, for simplicity, the fate of 
a coherent radial oscillation in the true minimum <p = (/„., 0) of the potential 
(12). Writing fa 1=1 /w + X an d 4> — the angular (it) oscillations are governed 
by the equations 113 



with m 2 = 2Xf 2 and A n = 4(mJ + fc 2 )/m^, gv = 2xo//tt- The Lame equation 
(13) has unstable, i.e. exponentially growing solutions in a wide range of 
parameters (A w , q n ). The unstable regions focus near A„ — n 2 (n = 1, 2, . . .) 
when q-ir — > 0, but broaden as q„ grows. 

The quantum field theoretic description of this process has recently been 
worked outt£l It turns out that the final state is described as a squeezed state 
of coherent pion pair excitations in the spontaneously broken chiral vacuum. 
The classical equation (13) governs the amplitude functions determining the 
squeezed state. The numerical results show a strong peak of the pion number 
in the momentum range where the Lame equation (13) exhibits instability. 
Strong charge and momentum correlations among pions are also found. For 
pions of opposite momentum these correlations are dictated by those arising 
from the decay a — > tttt while those among pions with nearly equal momenta 
are predominantly given by quantum correlations among identical particles 
(HBT correlations). 

4 The Anomaly Effect 

In the presence of strong electromagnetic fields, the Adler-Bell-Jackiw anomaly 
describes an additional interaction that affects only neutral pions. The inter- 



—2 + /Lr - 2q ff cos(2z) n k (z) = 



(13) 



where the radial (a) oscillations are approximated by 



X{t) = XoCOs(m CT i + a) 



(14) 
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action has the form 

5C = — %-{E-3)fa (15) 
shifting the minimum to an oblique position on the chiral circle given by 

(foUin = ~^^E ■ H. (16) 

This shift provides a bias in the dynamics of the pion field in the ir° direction 
whose sign depends on the sign of the invariant E ■ H. For heavy ion collisions 
at a nonvanishing impact parameter b ^ 0, this-sign differs between the regions 
of space above and below the scattering planell3 

The effect of the anomaly is numerically small even in strong electromag- 
netic fields due to the presence of the factor a in (15). However, its effect is 
not negligible because of the spatial homogeneity of the interaction (15) over 
distances r 3> m~ . Consider a highly energetic collision along the z-axis of 
two equal, heavy nuclei at an impact parameter b = Ra, where Ra is the 
nuclear radius. The pseudoscalar electromagnetic invariant is then given by 

- - 2Z 2 e 2 -y 2 

E - H =jwwm x " {bxvl (17) 

where R\i 2 = 1 2 { Z ± vt) 2 + {x±_ ± \b) 2 . At high energy (7 ^> 1) the value 

of E ■ H is strongly peaked near t — 0, and the anomaly interaction can be 
approximated by a quasi-instantaneous kick to the ir° field at the moment of 
encounter of the two nuclei. As noted above, the kick is small but spatially 
homogeneous over a region of the size of the nuclear radius Ra- 

To investigate the effect of the anomaly kick-cp the formation of DCC 
domains, we have performed numerical simulationsEZI of the linear sigma model 
in the modified "quench" scenario Ell with a schematic initial kick 

A<j> 3 = sgn(y)a n ml, (18) 

where y = denotes the scattering plane. According to (17) the value of a„ 
should lie in the range 0.05 - 0.1, for a value of 7 = 100 corresponding to 
nuclear collisions at RHIC. Our simulations show that the anomaly kick has 
a surprisingly large effect on the time evolution of the ir° field, although the 
initial kick is tiny when compared with the fluctuations of the pion field on 
length scales of order m^ 1 . The effect and the asymmetry with respect to 
the scattering plane is best revealed by looking at the time-like component of 
the neutral axial current, Aq, as a function of time (see Fig. [j]). While the 
sign of Aq varies over time, the sign always differs between the upper and 
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lower half-space. We also found that no other component of the axial or the 
vector current shows a similar behavior, and we have confirmed that the effect 
remains unaffected by a restriction of the anomaly kick to a finite region of 
space. 

ft appears difficult to detect this polarization of the axial charge directly, 
because pionic observables are usually not sensitive to the sign of Aq. However, 
we have found that the effect is sufficiently strong to slightly distort the isospin 
symmetry of DCC formation. This effect could be large enough so that it may 
be observable in a large sample of events, by comparing the average number 
of produced neutral pions to the average number of charged pions (see Fig. ||). 
The enhancement is concentrated at small momenta (k < m n ) and hence may 
be extracted by averaging the tt° /ir^ ratio over many events for high and low 
pion momenta, respectively. 
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Figure 1: <Aj) up 



(^(^uppcr, and 
(^o)lowcr as a function of proper time. The 
calculation is for an anomaly strength a n = 
0.1 and trip-average is taken over 10 events. 
[From rcf.EZf 



Figure 2: The solid line is the proper time 
evolution of N(ir°)/N(n) averaged over 10 
events. The dashed line is the same except 
that the momentum integration is limited 
to fc < 250 MeV. [From ref.EJ] 



5 Summary 

Coherent excitations of the chiral order parameter field, i.e. the physical pion 
field, in high energy heavy ion collisions can generate unusual distributions of 
the charge ratios of produced pions as well as peculiar pion charge correlations. 
The three-dimensional expansive cooling of an equilibrated fireball with initial 
temperature T > T c , the critical temperature of the chiral phase transition, 
provides a "natural" and plausible mechanism of such excitations in the form 
of domains of disoriented chiral condensate (DCC). Since the condition of 
approximately vanishing isospin is crucial for the phenomena discussed here, 
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they are most likely to occur in the nearly baryon-free central rapidity region 
in nuclear reactions at RHIC. 

The electromagnetic axial anomaly enhances the spectral power of the ir°- 
excitations at small momentum and manifestly breaks the tt°/it ± symmetry. 
This effect may result in a neutral pion ratio N(ir°)/N(ir) > -| in this region 
of phase space. The anomaly effect could also lead to an overall enhancement 
of the formation of extended DCC domains in semiperiphcral collisions of very 
heavy nuclei at RHIC or LHC. 
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